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ABSTRACT 

A study and pursuant development is described of a digital computer 
program for the computation of transfer function sensitivities in a 
symbolic form suitable for storage and subsequent repetitive numerical 
evaluation over a range of frequencies or parameter values . The algo- 
rithm is implemented in conjunction with the Network Analysis for Systems 
Applications Program (NASA?) version developed by R. S. Schwartz at 
Northeastern University, Several example problems, as well as sugges- 
tions for possible improvement of the algorithm and its extrapolation 
into the areas of optimization and automatic design, are included. 
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I. INTRODUCTION 



A. STATEMENT OF THE PROBLEM 

The purpose of this research was twofold: (1) to investigate the 

feasibility of providing for the automatic digital machine computation 
of transfer function sensitivities in a symbolic form suitable for 
storage and subsequent numerical evaluation over a range of frequencies 
and parameter values or for use by a minimization/optimization procedure 
as a readily available source of gradient data, and (2) to implement 
such a program, if feasible, in conjunction with the Network Analysis for 
Systems Applications Program (NASAP) . 

There were formidable obstacles encountered in the course of the 
research which may best be illustrated by a preliminary examination of 
the computer core storage required for implementation of the proposed 
procedure if it were meshed with an existent circuit analysis program 
such as NASAP. The NASAP program is capable of manipulating up to 10,000 
independent loop products (generated through use of topological flowgraph 
techniques) , once as a numerator polynomial in ^ and once as a denominator 
polynomial; that is, by iterative use of a single vector variable, LOOP , 
dimensioned to 10,000. Each loop product may represent from one to 
twenty-three constituent elements, and the analytical expression for the 
partial derivative of a quotient involves four separate components - 
numerator polynomial, denominator polynomial, partial derivative of the 
numerator and partial derivative of the denominator. The implication, then, 
is a storage allocation requirement of four vector variables each of 
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dimension 230,000 words (920K bytes). This, of course, would provide 
for the sensitivity with respect to a single parameter; if the gradient 
is desired the requirement for the two partial derivative variables 
would be multiplied by a factor of 23 again, for a total of 11,500,000 
words. Consider now an even more staggering statistic. Since the NASAP 
program, like most programs written in the FORTRAN language, assigns 
only one alphanumeric character to a word (four bytes in the IBM 360/67) , 
there would be an additional factor of four; three for the element 
identifying characters and one for a code to indicate the algebraic 
formulation required. This initial examination, then, would indicate the 
need for an impossibly high storage allocation of 46,000,000 words 
(184,000K bytes!). Some degree of success was achieved in the reduction 
of this figure to a manageable size, but at the expense of limitations 
on the original objective. The extent of the progress and associated 
restrictions will become clear in later sections. 

B. DISCUSSION OF THE NASAP PROGRAM 

Prior to discussing the sensitivity algorithm in any detail, it 
would perhaps benefit the reader to summarize the characteristics and 
capabilities of NASAP since it is upon this foundation that the proposed 
technique is firmly based. 

The revised version of NASAP is a digital computer program for the 
flowgraph analysis of lumped, linear active electrical networks designed 
to provide any desired circuit transfer function in either symbolic or 
numerical form. The program utilizes the "mixed** method to generate 
circuit equations, i. e., operations are performed on a tree-branch, 
link-branch topology so that Kirchoff’s voltage law may be utilized to 
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express each link voltage in terms of tree voltages and Kirchoff’s 
current law may be used to express each tree current in terms of link 
currents. These relationships, in association with volt-ampere and source 
dependency relationships, comprise a formulation defining a canonical form 
called a primitive flowgraph [3]. 

The user of NASAP prepares a circuit for analysis by: (1) converting 

to a linear model if necessary, (2) labeling each circuit element by type 
(R,L,C,V,I) and a two digit number V7hich uniquely identifies it, and (3) 
numbering, each node as desired except for the ground node which must be 
designated node 1, The user then inputs this connectivity data in the 
following format: 

ELEMENT FROM TO ELEMENT CONTROLLING PARAMETER 

NAME NODE NODE VALUE (For controlled source) 

Following this are coded sets of data describing the transfer functions 

desired, in the form of the ratio of a response variable (voltage or 

current associated with some circuit element) to an excitation variable 

(voltage or current), and the forms of program output desired. The NASAP 

tree-link hierarchy [20] and the order in which the elements were input 

to the program determines which elements are placed in the tree, thus 

providing the user with some flexibility regarding its choice.- Given the 

topology, the machine then forms the fundamental outset matrix and from 

this constructs the primitive flowgraph. 

A transfer function between any two nodes of the primitive flowgraph 
may be obtained by adding a "closure” [20] and solving the resulting 
"closed" graph by application of Mason’s gain formula in the following 
manner. If the transfer function between node j and node k of a flow- 
graph is equal to Tjj^ and if the graph is then closed by the addition of 
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a branch from k to j whose transmission gain is equal to the reciprocal 
of then the loop gain with the added branch is equal to unity and 

the determinant of the flowgraph vanishes. This reduces the solution of 
T to the problem of summing loop-set transmittance products, a procedure 

highly suitable to machine solution. Readers desiring more detailed 
information concerning this procedure are directed to Ref. 20. 

The procedures utilized by NASAP for labeling, sorting and storing 
each loop-set consist primarily of: 

(1) As each loop-set is determined by the search routine, a signed 
integer label is generated from which the loop-set may be 
uniquely regenerated . 

(2) Loop-sets with the same label but opposite sign are eliminated 
to reduce storage requirements. 

(3) Loop-sets are sotted and stored in order of decreasing magnitude. 

(4) The completed loop-set array represents the symbolic solution of 
either the numerator or denominator polynomial in coded form. 

From successive solutions the transfer function may be evaluated 
either numerically or symbolically by decoding and summing each 
loop. 

C. A COMPARISON OF TWO FAVORABLE TOPOLOGICAL METHODS 

A topological method other than the flowgraph procedure just described 
is the so-called ”k-tree" approach [6] which is based on Maxwell's topolo- 
gical formulas for node voltage equations. That is, transfer functions are 
derived from the evaluation of determinants formed from circuit nodal 
equations ;however , they may equivalently be formed in a topological fashion 
by summing all tree admittance products and appropriate 2-tree, 3-tree, ... 
k-tree admittance products. 
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In both methods, k-tree and loop-set, the solution of passive 
circuits is obtained as a sum of linear products of the network parameters, 
thus eliminating the requirement for operations with polynomials, and 
facilitating machine evaluation of symbolic solutions. In addition, the 
two methods compare favorably with respect to the number of loop-trans- 
mittance products and tree-admittance products for a given circuit, since 
the maximum number of trees possible for a given topology is only one 
greater than the maximum number of loop-sets contained in the corresponding 
primitive flowgraph [20]. Each tree-admittance product is unique; however, 
in the flowgraph method several loop-sets are frequently found to have 
identical transmittance product magnitudes, although only one will remain 
after cancellation. Thus each remaining product is unique. A primary 
advantage of both procedures is this feature which obviates the evaluation 
of zero terms, as they would be if the system determinant were evaluated 
by matrix algebra. 

Although the two methods appear equal from the above discussion, a 
very important difference arises in the solution of active networks since 
the k-tree approach requires the graph to be modified and divided into 
separate current and voltage graphs from each of which the complete k-trees 
must determined, whereas in the loop-set analysis the procedure is invariant 
for passive or active networks. 

Another significant disadvantage of the k-tree analysis is its 
dependence upon the particular transfer function, to be evaluated, i.e., 
impedance, admittance, gain-, etc., each requiring different types of k-tree 
searches. Thus, the number of k-tree searches will vary with the number of 
transfer functions as well as the type desired. Therefore, the flowgraph 
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method appears to have a distinct advantage in this respect, since the 
procedure remains identical regardless of the number or types of transfer 
functions. For these reasons the flowgraph approach was chosen for imple- 
mentation in NASAP, 

D. SUMMARY 

It was the objective of this study to determine the feasibility of 
utilizing the loop-set symbolic form of the transfer fmnction originated 
in the manner described in (B) above as a basis for computing and storing 
sensitivity information in similar symbolic form, thereby vastly simplify- 
ing the repetitive computation of accurate network sensitivities and 
rendering such solutions free of the difficulties involved with differen- 
tial, or numerical techniques. 
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II. DEVELOPMENTAL PROCEDURE 



Utilizing the revised version of NASAP developed for the PDP-10 and 
CDC-3300 computers by R. S. Schwartz [20], adaptation to the IBM 360/67 
computer at the Naval Postgraduate School, Monterey was accomplished by: 

(1) altering program statements necessary to reconcile the difference in 
word lengths between the two types of computers, (2) making various 
syntax changes necessitated by technical differences in the implementa- 
tion of FORTRAN on the different machines, and (3) changing the modular 
form and adding program statements necessary to provide for the execution 
of the algorithm as intended by the author. Details of the conversion are 
apparent from comparison of the program listing available at the W. R. 
Church Computer Center, USNPGS, Monterey, California with that of Ref. 20. 
The listing found in the appendix of this thesis contains only the NASI, 
NASA, and SENS subroutines. The first two are included since extensive 

A 

modifications were incorporated in the course of implementing the 
sensitivity feature, and the latter because it is entirely new. 

As previously described, the NASAP program has available, in symbolic 
form, the uncancelled loop products created from a topological flowgraph 
analysis of a given circuit. It was from this base that the writer 
devised a program to compute circuit transfer function sensitivities in 
an algorithm developed as follows. 

Since, in subroutine NASA of the original program, there implicitly 
exists a requested transfer function, 

T(s) = K HifI in which, 

D(s) 
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N(s) = bgs’^ + + . . . + b^sP + . . . + bj^pj^sq 

D(s) = a„s’^ + a,s™ + . . . + a.s^ + . . . + a 

0 1 1 LPD 

where n, m, p, q, p', q’ , ... are integer exponents of the Laplace 

variable not sorted into any ascending or descending order and with 
repetitions possible; and where a^, b^ represent the i^^ loop products 
(algebraic product/quotients of circuit parameters in symbolic form); 
and where LPD is the number of uncancelled loop products in the denomi- 
nator polynomial of T(s) and LPN the number in the numerator. Each loop 
product is actually determined from one 32-bit computer word, wherein 
the presence (1) or absence (0) of each element is indicated by the 
corresponding bit value. This, of course limits the number of circuit 
parameters which can be accommodated by the program to NBITS-1 = 31, 
although in the current version this upper limit is set to 23 due to other 
considerations , 

It was desired, then, to compute the sensitivity of T(s) in symbolic 
form, to a specified circuit parameter P^^, Since sensitivity is essen- 
tially a rate-of-change, it follows that the transfer function sensitivity 



may be defined as: 



DCs).i^-N(s).^ 



8P, 



= -9T(s) = 8(N(s)/D(s)) _ 

k k k 



D^(s) 



( 1 ) 



To implement this fundamental approach utilizing the data available 
in NAS4, the following simple algorithm was developed, (See Flow-Chart, 
Fig. 1): 

(I) Observe the j factor of the i loop product, b. . , in the numerator 



of T. 



th 



Is bj^j in the numerator of the i loop product, b^? 
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YES: 



(a) Is this factor the parameter 



YES: (i) Set flag to indicate that b. contains P, , 

1 K. 



(ii) Form 



NO; Form 



(b) N. . = b. , 

13 13 



V k^ij 
(^k/11 



NO: 



(a) Is this factor Pj^? 



YES : (i) Set flag 

(ii) Change sign of 






(iii) Form 



/3N \ 

' i^X' 



NO: Form 



(SI 



k/ij 

= 1/b. , 

13 



k' i 






(b) N. .= 1/b. . 

13 13 



(II) Repeat step (I) for the factor of the loop product. 

Continue until all factors are exhausted (j>NFN(i), where NFN(i) 

represents the number of factors in the i^^ loop product of the 
numerator polynomial) . 

(Ill) Does j contain Pj^? (is flag s6ttol?) 

' k 'i 

YES: Continue 



NO: Set 



:f^\ = 0 

hPv/ ,• 



IV) 



k" i 

(A) i = i + 1 

(B) Repeat (I) (II), and (III) for the (i+l)^“ loop product. 
Continue until all loops are exhausted (i>LPN) . 



xth 
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(V) Repeat the entire procedure beginning with step (I) for the 
denominator of T(s). 

At this point the sensitivity is implicitly available in the desired form 
(eqn. 1), where: 



D = 


LPD 

Z 

i=l 


NFD(i) 

n 

3=1 


EXPDS(i) 

D. . s 
13 


N = 


LPN 

i=l 


NFN(i) 


EXPNS(i) 

N.js 


8D 


LPD NFD(i) 

= E n 
i=l j=l 


(^Ji3 


9N 


LPN 

Z 

i=l 


NFN(i) 

n 

j=i 


/8N '\ 



NFN(i) = Number of factors in the i^^ loop product of numerator 
polynomial 

NFD(i) = Number of factors in the i^^ loop product of denominator 
EXPNS(i) = Exponent of ^ in the i^^ loop product of numerator 
EXPDS(i) = Exponent of ^ in the i^^ loop product of denominator 
One feature of the programming which became necessary in order to 
reduce storage requirements for the symbolic sensitivity component variables 
SNUM, SDEN, PARTLN, and PARTED was a packing scheme which in effect 
reduced vector sizes by a factor of one-fourth. In the NASAP program, as 
in most FORTRAN programs, each alphanumeric character is represented by 
a four-byte word, wasting a great deal of core storage since one byte is 
sufficient for one character. This exemplifies one of the most unattractive 
aspects of the FORTRAN language, especially when dealing with character 
strings or literals as in the case of the NASAP symbolic version. 

A remedy for the problem which was utilized in this program is a 
very short but efficient packing subroutine which has four arguments: 
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(1) a four-byte variable to be “packed** 

(2) an integer constant specifying the byte number of the variable 
in (1) to be assigned the value of the packing variable 

(3) packing variable 

(4) an integer constant specifying the byte number of the variable 
in (3) which is to be assigned to the variable in (1) 

Upon calling the pack subroutine, dummy variables, locally declared 
LOGICAL*! (one byte) , receive data passed from the specified byte of the 
packing variable, assign this to the specified byte of the variable to be 
packed, and then return this new **packed** version to the calling program. 

In this manner the sensitivity component variables were packed with para- 
meter name (e.g. LOl, R04, etc.) in bytes one through three, and byte 
number four was packed with a code number to indicate the form in which 
the parameter should appear. This code appears below: 

0 — Numerator or denominator = z, partial = z 

3 — Numerator or denominator = z, partial = 1 

4 — Numerator or denominator = z, partial = 0 

6 — Numerator or denominator = 1/z, partial = 1/z 

7 — Numerator or denominator = 1/z, partial = -1/z^ 

9 — Numerator or denominator = 1/z, partial = 0 

With this scheme the program was rendered capable of processing net- 
works of up to 200 independent loops to provide for ‘sensitivity with 
respect to a single parameter, or by successive solutions to provide for 
any or all sensitivities desired. With this loop-set limit, the entire 
program, including NASA? requires 180K of core. A guide for core storage 
allocation may be found in the appendix for circuits of greater complexity. 
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Subsequent to the packing operation, data is transferred to a new 
subroutine SENS for manipulation into the desired form and for numerical 
evaluation. In this operation it was found impractical, of course, to 
permit each product element to occupy a separate storage location in the 
computer, since this would essentially require dimensions on the order 
of the square of the maximum number of loop products. Therefore, numerical 
evaluation of the constituent components of the sensitivity are recursively 
computed and simultaneously combined at this point in the program to form 
the numerator and denominator arrays, PARTI and PART2. Eliminating all 
zero vectors and combining like powers of ^ reduced storage requirements 
for this operation essentially to twice the order of the denominator 
polynomial. 
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III. USER'S RULES 



Users of the sensitivity subroutine must be familiar with the simple 
NASAP coding rules contained in Refs. 1 and 2 as modified below: 

(1) Only one *'OUTPUT" and/or "ROOTS** transfer function request should 
be made in a given problem since the sensitivity subroutine is 
designed to store data relevant to only one numerator polynomial 
at a time. Of course, the user may wish to follow one set of 
problem cards with another. 

(2) The **SYMBOLIC SOLUTION** card must be inserted following **EXECUTE**. 

(3) No **MODIFY** requests are permitted. 

(4) The card following **SYMBOLIC SOLUTION** must contain the words 
(only the underlined portions are essential for machine 
recognition) ; 

SENSITIVITY PARAMETE R=L02 (or ROl, C09 , etc.) 

Any number of blanks or identifying characters may be included 
so long as the underlined portions appear somewhere in columns 
1-80 in the order shown. Notice that the parameter (with respect 
to which the sensitivity is desired) is indicated by the consecu- 
tive characters, (=) , (R, L, or C) , ( two digit element number). 

(5) If it is desired to have the sensitivity magnitude and phase 
evaluated, the next card in sequence must contain a maximum 
frequency in columns 1-11 and a minimum frequency (not zero) in 
columns 21-31, preferably both in standard FORTRAN formal Ell. 4, 
although normal floating-point form is satisfactory ^ all 
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significant digits (plus sign) remain XN^ithin the specified 
11 columns. Specifying these frequencies limits v/ill result in a 
print-out of the sensitivity magnitude and phase at 200 evenly 
distributed frequencies in the interval. 

( Caution ; If the user omits this option and desires to follovzs 
the current problem deck with others, a blank card must be 
inserted in place of the frequency interval card) . 

An interested user may provide for the plotting of the generated 
magnitude and phase data by changing the variables RMAGN, ANGLE, and 
FREQ into array variables in the SENS subroutine, then passing them to 
an existing library plotting subroutine. This will, of course, consider- 
ably increase the core storage required for the program. 
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IV. INTERACTIVE APPLICATIONS 



The NASAP program, in conjunction with the sensitivity subroutine, 
may be easily adapted for use under a time-sharing system such as CP/CMS 
with only a few simple statement changes which depend upon the user’s 
choice of file structure. One choice successfully used at the Naval 
Postgraduate School is to remove the four initial statements: (1) CALL 

NASAP, (2) STOP, (3) END, and (4) SUBROUTINE NASAP. It is possible, then, 
to "offline read" the program into CP in the following modular form: 



FILE 


PROGRAM CONTENT 


RSS 


Program RSS (Driving program) 


NASI 


NASI, FUNCTION ITRAI^I 


NAS 2 


NAS 2 


NASS 


NASS 


NASA 


NASA, SUBROUTINE PACK 


NASS 


NASS 


SENS 


SENS, BLOCK DATA 



Under CP/CMS, the entire program can then be implemented by executing 
"RSS", with data being input at the terminal, or if the three "READ" 
statements in NASI are changed appropriately, data may be read from a 
separate file so defined. 
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V. DISCUSSION AND RECOMMENDATIONS FOR FURTHER IMPROVEMENT 



As the reader may have gathered, the feasibility of computing sensi- 
tivities of a transfer function to all circuit parameters and storing 
them in symbolic form appears doubtful under present storage limitations. 
What follows, however, is a recommendation for further reducing storage 
requirements somewhat while maintaining the symbolic concept. If only 
a single numerical sensitivity were desired, it would of course be 
considerably more efficient to recursively evaluate the sensitivity 
from successive loop-sets with their appropriate values first as given, 
then in partial derivative form. This procedure would eliminate most 
of the storage problems. The singular advantage of the symbolic approach, 
however, is that evaluation of the sensitivity function at various fre- 
quencies, of for various parameters, is a very simple matter of substi- 
tution, rather than the use of time-consuming iterative or numerical 
techniques which may pose stability or accuracy difficulties. 

It is apparent that substantial changes are indicated in the proposed 
algorithm (or perhaps only the techniques for implementing the algorithm 
which will permit the sensitivity subroutine to cope with circuits of 
greater complexity by means of a reduction in the large storage require- 
ments. However, it must be remembered that the resulting exponential 
increase in the number of independent loops with respect to a concomi- 
tant increase in the number of nodes [20] will in itself limit the 
complexity of topological/f lowgraph methods, particularly in maximal- 
planar networks. Indeed, the difficulties of dimensionality plague all 
circuit analysis programs. 
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One suggested method of reducing storage requirements would be 
to extend the previously mentioned concept of utilizing individual 
bits of the computer word to represent constituent elements of the 
net\^ork. Since in this case it is necessary to know which of six 
mathematical formulations is appropriate for each element of the 
numerator and denominator polynomials, instead of only two formula- 
tions as in the case of transfer function loops, it would appear 
that a scheme requiring only three variables (2^ = 8 possibilities), for 
a numerator and denominator each dimensioned to the maximum number of 
loops (10,000 in the present version) would suffice to allow sensitivi- 
ties to be computed and stored in symbolic form. This would repre- 
sent a substantial reduction from the 460,000 that are presently required 
for a 10,000-loop circuit to accomplish the same task. This suggestion, 
in effect, offers an execution time-vs-storage tradeoff, since a coding/ 
decoding algorithm must be developed which will intrinsically involve 
a tremendous number of manipulations. In addition it would require an 
almost total rewrite of the FORTRAN programming, particularly that 
portion meshed with the NASA? NASA subroutine, and it may very well 
lead to difficulties not apparent without further investigation. 
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APPENDIX A SAMPLE PROBLEMS 



(I) Sample problem No. 1 - Emitter Follower 



(a) Emitter follower linear model: 




(b) NASAP coding: 



NASAP PROBLEM NO, 1 

♦EMITTER FOLLOWER 

VI 1 2 1 

1114 1 

R1 2 3 100 

R2 3 4 IM 

R3 3 1 1.5K 

R4 4 1 100 

Cl 3 4 6PF 

C2 3 1‘200PF 

12 4 1 0.04VC2 

OUTPUT 

VR4/ VVl 

EXECUTE 

SYMBOLIC SOLUTION 
SENSITIVITY PARAMETER=R04 
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(c) Sensitivity output for emitter follower: 



numerator POLYfJOMIAL : SENSITIVITY OF VR04/VV01 TO 
0.6944E 40 TIMES S** 0 

0.1042E 31 TIMES S** -1 
0.6667E 17 TIMES S** -2 



DENOMINATOR 


0. 2778E 


19 


0. 3473E 


31 


0. 1085E 


43 


O.llllE 


30 


0.2845E 


16 



POLYNOMIAL : 
TIMES S** 2 

TIMES S** 1 

TIMES S** 0 

TIMES S** -1 
TIMES S** -2 



SENSITIVITY 



R04 
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(II) Sample Problem No. 2 - Butterworth Filter 



(a) Original linear circuit: 




(b) NASAP coding: 



NASAP PROBLEM 

♦THIRD ORDER BUTTERWORTH 

R1 4 2 1 

R2 3 1 1 

Cl 2 1 1 

LI 2 3 2 

C2 3 1 1 

VI 1 4 1 

OUTPUT 

VR2/VV1 

EXECUTE 

SYMBOLIC SOLUTION 
SENSITIVITY PARAMETER=R02 
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(c) Sensitivity output for Sample Problem No. 2 - Butterworth 
Filter: 



NUMERATOR POLYNOMIAL : SENSITIVITY OF VR02/VV01 TO 
0.5000E 00 TIMES S** 0 

0.5000E 00 TIMES S** -1 

0. 2500E 00 TIMES S** -2 



DENOMINATOR POLYNOMIAL ; SENSITIVITY 
O.IOOOE 01 TIMES S** 4 

0.4000E 01 TIMES S** 3 

0.8000E 01 TIMES S** 2 

O.IOOOE 02 TIMES S** 1 

0.8000E 01 TIMES S** 0 

0.4000E 01 TIMES S** -1 
O.IOOOE 01 TIMES S** -2 



R02 
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(Ill) Sample Problem No. 3 - Transformer circuit 
(a) Original circuit with ideal transformer: 




(b) Linear model: 



© 




(c) NASAP coding: 



nasap problem 

♦INPUT IMPEDANCE 

11 1 2 1 
R1 2 4 3 
Cl 3 1 0.5 
LI 2 3 0.25 

12 3 1 -2 IRl 
V2 4 3 -2 VCl 
OUTPUT REQUEST 
VI 1/ I 1 1 
EXECUTE 

SYMBOLIC SOLUTION 
SENSITIVITY PARAMETER=L01 

.1592E 03 .1592E 00 
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(d) Sensitivity output for Sample problem No* 3 "" Transformer 



circuit : 



NUMERATOR POLYNOMIAL : 



0. lOOOE 


01 


Tl 


IMES 


S** 


-0.3333E 


00 


TIMES 


s** 


0.2778E 


-01 


Tl 


IMES 


s** 



DENOMINATOR 


POLYNOMIAL 


0.1778E 


01 


TIMES 


s** 


0. 2667E 


01 


TIMES 


s * * 


0. 1444E 


01 


TIMES 


s** 


0.3333E 


00 


TIMES 


s * * 


0.2778E- 


01 


TIMES 


s * * 



ITI VI TY OF VI 01/ I I 01 TO LOl 



ENSITI VI TY 



SENS 

1 

0 

1 

: S 

2 

1 

0 

1 

2 
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(_e) Sensitivity magnitude and phase output for 



No. 3 - Transformer circuit: 



evaluation of sensitivity 



frequency 

0.1592E 00 
0.9584E 00 
0.1758E 01 
0.2557E 01 
0.3356E 01 
0.4155E 01 
0.4954E 01 
■0.5754E 01 
0.6553E 01 
0.7352E 01 
0.8151E 01 
'0.8950E 01 



MAGNITUDE 

0.4352E 00 
0.9268E-01 
0.5082E-01 
0.3498E-01 
0. 2666E-01 
0.2154E-01 
0.1806E-01 
0.1556E-01 
0. 1366E-01 
0.1218E-01 
0. 1098E-01 
O.lOOOE-01 



U.b23UE 02 
0.6330E 02 
0.6409E 02 
0.6489E 02 
0.6569E 02 
0.6649E 02 
0.6729E 02 
0.6809E 02 
0.6889E 02 
0.6969E 02 



0.1432E-02 
0. 1414E-02 
0.1397E-02 
0.1380E-02 
0.1363E-02 
0.1346E-02 
0.1330E-02 
0.1315E-02 
0. 1300E-02 
0.1285E-02 



0.1480E 03 
0.1488E 03 
0.1496E 03 
0.1504E 03 
0.1512E 03 
0.1520E 03 
0.1528E 03 
0.1536E 03 
0.1544E 03 
0.1552E 03 
0.1560E 03 
0.1568E 03 
0.1576E 03 
0.1584E 03 
0.1592E 03 



0.6049E-03 

0.6016E-03 

0.5984E-03 

0.5952E-03 

0.5921E-03 

0.5890E-03 

0.5859E-03 

0.5828E-03 

0.5798E-03 

0.5768E-03 

0.5739E-03 

0.5709E-03 

0.5680E-03 

0.5652E-03 

0.5623E-03 



Sample Problem 



PHASECIN RADIANS) 

0.1763E 00 
-0.1267E 01 
-0.1405E 01 
-0.1457E 01 
-0.1484E 01 
-0.1501E 01 
-0.1512E 01 
-0.1520E 01 
-0.1526E 01 
-0.1531E 01 
-0.1535E 01 
-0.1538E 01 



-0.1566E 01 
-0.1566E 01 
-0.1566E 01 
-0.1566E 01 
-0.1566E 01 
-0'. 1566E 01 
-0.1566E 01 
-0.1567E 01 
-0.1567E 01 
-0.1567E 01 



-0.1569E 01 
-0.1569E 01 
-0.1569E 01 
-0.1569E 01 
-0.1569E 01 
-0.1569E 01 
-0.1569E 01 
-0.1569E 01 
-0.1569E 01 
-0.1569E 01 
-0.1569E 01 
-0.1569E 01 
-0.1569E 01 
-0.1569E 01 
-0.1569E 01 
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(IV) Sample Problem No. 4 - I.F. Transformer 



(a) I.F. Transformer circuit: 




(b) Linear model: 




(c) NASAP coding: 



NASAP PROBLEM 
*I.F. TRANSFORMER 
II 1 2 1 
R1 2 1 .5M 
Cl 2 1 245PF 

R2 2 3 10 
LI 5 4 .49MH 
L2 4 1 .OIMH 
L3 4 5 .49MH 
R3 5 6 10 
C2 6 1 240PF 
C3 2 6 2 PF 
R4 6 1 .5M 
OUTPUT 
VR4/ I I 1 
EXECUTE 

SYMBOLIC SOLUTION- 
SENSITIVITY PARAMETER=C03 
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(d) Sensitivity output for Sample Problem No. 4 - I.F. 



Transformer : 



numerator polynomial : 



0.4720E 


18 


Tl 


IMES 


S** 


0.4786E 


16 


Tl 


IMES 


s** 


0.1293E 


14 


Tl 


IMES 


s** 


0 . 4049E 


10 


Tl 


IMES 


s** 


0.4725E 


06 


Tl 


IMES 


s** 


0.2309E 


02 


Tl 


IMES 


s** 


0.4000E- 


•03 


Tl 


IMES 


s** 



DENOMINATOR 


POLYNOMIAL 


0.4760E 


42 


Tl 


IMES 


s** 


'0.4827E 


40 


Tl 


IMES 


s** 


0.1304E 


38 


TIMES 


s** 


0.4095E 


34 


Tl 


IMES 


s * * 


0.4805E 


30 


Tl 


IMES 


s * * 


0.2375E 


26 


Tl 


IMES 


s** 


0.4266E 


21 


Tl 


IMES 


s * * 


0.4074E 


15 


TIMES 


s * * 


O.IOOOE 


09 


Tl 


IMES 


s** 



ITIVITY OF VR04/I I 01 TO C03 



ENSITI VITY 



SENS 

5 

4 

3 

2 

1 

0 

1 

: S 

6 

5 

4 

3 

2 

1 

0 

1 

2 
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(V) Sample Problem No. 5 - Wien-Bridge circuit 



(a) Original circuit: 




(b) Linear model: 




(c) NASAP coding: 



NASAP PROBLEn 
*V/I EN-RRIDGE CIRCUIT 
1112 1 
R1 2 1 1 
Cl 2 1 1 
R2 2 3 1 
C2 5 4 1 
V2 1 4 1 VRl 
OUTPUT 
VV2/ I I 1 
EXECUTE 

SYMBOLIC SOLUTION 
SENSITIVITY PARAMETER=R02 
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(d) Sensitivity output for Sample Problem No. 5 - Wien-Bridge 



circuit: 



NUMERATOR POLYNOMIAL : SENSITIVITY OF VV02/II01 TO 
IDENTICALLY ZERO 



DENOMINATOR 

O.IOOOE 01 

0.4000E 01 
0.6000E 01 

0.4000E 01 

O.IOOOE 01 



POLYNOMIAL : 
TIMES S** 4 

TIMES S** 3 

TIMES S** 2 

TIMES S** 1 

TIMES S** 0 



SENSITIVITY 



R02 



34 






1 - 
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I 

I 
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(VI) Sample Problem No. 6 - High-Q Filter 



(a) Filter circuit: 




(b) NASAP Coding: 



NASAP PROBLEM 
*HIGH - Q FILTER 
VI 1 2 1 
R1 2 4 20K 
R2 5 3 lOK 
R3 6 4 20K 
R4 6 7 100 
R5 7 5 100 
R6 5 1 2K 
R7 7 1 IM 
Cl 2 3 .02UF 
C2 5 4 ,04UF 
C3 6 3 ,02UF 
II 1 7 150 IR4 
1215 150 IR5 
OUTPUT 
VR6/VV1 
EXECUTE 

SYMBOLIC SOLUTION 
SENSITIVITY PARAMETER=C02 
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(c) Sensitivity output for Sample Problem No. 6 - High-Q filter: 



NUMERATOR POLYNOMIAL 
-0.2684E 09 TIMES S** 
0.3125E 12 TIMES S** 
0.3125E 09 TIMES S** 
0.1031E 06 TIMES S** 
.0.1250E 02 TIMES S** 
.0.5000E-03 TIMES S** 



SENSITIVITY OF VR06/VV01 TO 
6 
5 
4 
3 
2 
1 



DENOMINATOR 


0.1650E 


24 


0. 2869E 


21 


0.1887E 


18 


0.5643E 


14 


0.6908E 


10 


0.1575E 


06 


O.IOOOE 


01 



POLYNOMIAL : 
TIMES S** 6 

TIMES S** 5 

TIMES S** 4 

TIMES S** 3 

TIMES S** 2 

TIMES S** 1 

TIMES S** 0 



SENSITIVITY 



C02 
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APPENDIX B GUIDE FOR THE ALLOCATION OF CORE STORAGE 



* 






'k'k'k 



NUMBER OF LOOP-SETS POSSIBLE 
(VARIABLE DIMENSION) 




PROGRAM CORE PvEQUIRED 


* 200 




180K 


300 




205K 


400 




230K 


500 




255K 


600 




280K 


700 




30 5K 


** 800 




330K 


900 




35 5K 


1000 




380K 


5^** 5000 




505K 


10000 




630K 


Corresponds approximately 


to 


a maximal planar network 


of 5 nodes 






Corresponds approximately 


to 


a maximal planar network 


of 6 nodes 






Corresponds approximately 


to 


a maximal planar network 



of 7 nodes 
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PROGRAM LISTING 



SUBROUTINE NASI 

INTEGER EOFCKF , SENSFL, PARAM 

DIMENSION INLST( 80) ,NDDE( 21) 

COMMON/ SEN/ PAR AM (3 ) , S ENS FL , FRMAX , FRM I N 

C0MM3N/DATA/IVJSRK 100) , I NSR2( 24 , 6 ) , I TSKl ( 1 5 ,4) ,ITSK2(3 
1,6), ITSK3 

5(15,6) ,ITSK4( 10,5) , ITS K5 (9,10), I RAY (19, 2) , ISYMB(48), 

6 ITASK6( 17,12),I0PRTN(3,2) 

COMMON IND,NASA, MTRAN( 3) ,MD,IJTAG(15,2),IPRIN,Z(42 
1 ) , LL, MODNM( 23) , VALMOD( 23 , 1 0) , I LOOP (23 ) , L I , L I C , M9 , ZZ ( 23 
2 ) , 

3NAM( 23 ,3) , ITAG (5,8 ) , NMBRI6 ) , 

4M,LC0NC( 10,40) , NAME (23,2) ,MVI, I ST{ 40) , INC (38,42 ) ,N1,NA 
+ ,M2 

DATA IB/* • / , I E/IHE/ 

FRMAX=0. 0 
FRMIN=0 .0 

66 READ! 5,101 ,EN0=998) INLST 
65 DO 1 J=l,32 

IF( INLST (J )-IB )300 ,1,300 

300 DO 2 1=1,3 

IF ( INLST (J )-IOPRTN ( I, 1) )2, 350,2 

350 IF(INLST( J+1)-I0PRTN(I,2) ) 2, 360, 2 
2 CONTINUE 

GOTO 3 

360 G0T0(4,301 ,351 ) , I 

1 CONTINUE 
GO TO 66 

351 GO TO ( 67, 303, 352, 68, 68), NASA 

352 IPRIN=2 
IND=4 

WRITE (6,119) 

WRITE (6,116)INLST 

READ(5,101, END=1357) INLST 

IF( INLST( 1)-I0PRTN(3,1) ) 1357,1351,1357 

1351 IF ( INLST ( 2)- IE ) 1357, 1352, 1357 

1352 SENSFL=1 

DO 1353 K=3,80 

I F( INLST (K )-ISYMB(46) ) 1353,1355,13 53 

1353 CONTINUE 

1355 KP1=<+1 

DO 1354 I=KP1, 80 
I F( INLST( I ) . EO . IB) GO TO 1354 

1356 PARAM( 1) = INLST( I ) 

DO 13S0 MM=1,10 

IF( INLST(I+1).E0.1SYMB(MM) ) GO TO 1361 

1360 CONTINUE 

1361 PARAM(2) =MM-1 
DO 1362 MM=1,10 

IF( INLST(I+2).E0.ISYMB(MM) ) GO TO 1363 

1362 CONTINUE 

1363 PARAM( 3) =MM-1 

READ(5, 100, END=60) FRMAX, FRMIN 
GO TO 60 

1354 CONTINUE 

1357 SENSFL=0 
GO TO 60 

3 K= J 

GOTO (67, 69,68,68,69) , NASA 
67 WRITE (6,119) 

WRITE (6,109) 

NAS A=4 

68 WRITE ( 6,116) I NLST 
GO TO 66 

301 GO TO (57,303,305, 68, 68), NASA 
303 WRITE (6,304) 
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305 



40 0 



4 



5 



6 

7 

69 



10 



11 

12 

13 

14 

15 

16 

17 

18 

308 

19 



GOTO 68 
K=J+1 
Kl = 80 
K2=6 

DO 400J=1,23 
IL00P( J )=0 
MODNM( J )=0 
DO 400 L2=1,1D 
VALMODC J,L2)=0. 

Ll=6 
L2 = l 
M9=l 
ISGN=1 
LD = 0 
IZ=0 
LI1=0 
IS = 1 

WRITE (6,119) 

GOTO 69 
NASA=2 
Kl=30 
K2 =0 

WRITE (6,119) 

WRITE (6,116)INLST 
DO 5 1=1,21 
Z( I ) =0. 

Z( 1+21 )=0. 

NODE (I )=0 
DO 5 J=l,42 
INC( I, J )=0 
MTRAN(1)=0 
MTRAN( 2)=0 
MTRANI3 )=0 
MC=0 
N = 0 

IPRIN=0 

M=2 

M0VFL=23 
ISGN=+1 
LD = 0 
Ll=l 
M=M+1 

IF (M-MOVFL-3) 7,96,96 

L2=l 

IZ=0 

GO TO 66 

WRITE (6,116)INLST 

DO 47 J=K,K1 

ILST = ITRAN( INLSTI J) ) 

IFI INSRK ILST) )12,98,10 
L3 = INSR1 ( ILST) 

ILST=L3 

IF(INSR2(L3,L1))99,98,11 
L3=INSR2(L3,L1 )+l 
GO TO 13 
L3=l 

ILST=-INSR1 (ILST)-l 

GOTOI 14, 15, 16, 17,18,308) , LI 

L3=ITSK1(L3,L2 ) 

1=1 

GO TO 19 
L3=ITSK2(L3,L2 ) 

GO TO 19 
L3=ITSK3(L3,L2 ) 

GO TO 19 
L3=ITSK4(L3,L2) 

1 = 2 

GO TO 19 

L3=I TSK5(L3,L2 ) 

GOTO 19 

L3=ITASK6(L3,L2) 

GO TO (20,21,21,24,26,27,28,29,30,31,33, 
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+36,37, 38 » 39, 40 t 4 5, 46, 47, 51 , 6 6, 97, 98 , 99 ,41 , 42 ,43 , 44 , 
+310, 311, 31 2, 31 3, 31 4, 316, 316, 3 19, 316, 3 16, 316, 3 16, 321, 
+375) ,L3 

20 MN=10»MN+ILST 

21 I NC( MN+1 ,M) =4’!' ( L2-3 )*>!‘3/ ( ( L2-3 + 

IF(MN-N) 23,23,22 

22 N=MN 

23 NODE(MN) =NODE( MN)+1 
G0T0(46, 45,39) ,L3 

24 IF( ILST )99,47,25 

25 MN=ILST 
GO TO 46 

26 ISGN=-1 
GO TO 46 

27 INC( I,M) = 100’!=ISGN*ILST 
GO TO 46 

28 INCd ,M)=INC(I ,M) + ISGN^^MN 
GO TO ( 39,46), I 

29 INCH , M) =INC( I ,M) + ISGN-(10’!^MN+ILST ) 

GO TO 46 

30 LD=LD+1 

31 IZ=10=i'I Z+I SGN*ILST 
IF( ILST )99,47, 32 

32 L2 = L2+(120 + (L2-K2)=«=(2-{L2-K2)*6) )/36 

GOTO 47 

33 IF{ INC( 1,M)/ 100-4) 34,35,34 

34 LD=LD + 15 + ( ILST -10)* ( 5 82- ( I LST- 1 0 ) * ( 40+( ILST- 10) “S'! 22) ) ) 
1/100 

35 LD=LD+3-ILST 
IF(L1-6)322,315,99 

322 L2=6 

GOTO 47 

36 Ll=4 
L2=l 

ISGN=(4-ILST)/3 

I F( ( (I NC(1 ,M)/100)-4)/3 )46 ,98,46 

37 ISGN=( 4-IL ST)/3 
MC=MC+1 
MTRAN(MT)=MD 
MD=MC/2 
ME=MC-2*MD+1 
GOTO 46 

38 Ll=5 
M2=M-1 
MC=1 
MT=1 
GOTO 7 

39 L1=L1+1 
L2=l 
GOTO 47 

375 WRITE 16,376) 

NASA= 5 

40 JC=MC-2*MD+1 

IJTAG( MD,JC) =1 00*1 SGN*I LST 
GOTO 46 

41 JC=MC-2*MD+1 

I JTAG( MD, JC) =I JTAG ( MD, JC) + ISGN*MN 
GOTO (45, 46), ME 

42 JC=MC-2*MD+1 

I JTAGIMD,JC)=I JTAG(MD, JC)+ISGN*( 10*MN+ILST ) 

GOTO 46 

43 MT=2 

MTRANK 2)=MTRAM ( 1 ) 

GOTO 7 

44 MT=3 

MTRAN( 3)=MTRAM ( 1) 

GOTO 7 

310 IS=-2 
G0T0(45, 46) ,L2 

311 M0DNM(M9 )=100*ILST 
LI = 1 

L2=4 
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GOTO 47 

312 MODNMC M9) = M0DNM{ M9 ) +1 0+MN+ I LST 
GOTO 46 

313 M0DNH(M9)=M0DNM(M9)+MN 
GOTO 45 

319 IS=IS+1 

IF ( !S>!'LI+3 199,315,324 

314 IF(LI-10)315,315,326 

315 V=FLOAT(IZ) 

VALMODI M9,LI ) = V/10.0^«' (LD) 

L2=12 

IF( 15)331,99,332 

331 IF(LI-3)47,334,324 

334 IF(L3-36 >47,340,47 

340 IL00P(M9)=(VALM0D(M9,1 ) -V ALMO D( M9 , 2 ) ) /VALMODI M9, 3 )- 1. 0 

10001 

IF( IL00P(M9 ) ) 333,32 8,3 28 

332 IL00P(M9)=LI 

333 IF(LI1-IABS(IL00P( M9) ) ) 335,47,47 

335 LI1=IABS(IL00P (M9) ) 

GOTO 47 

316 IZ=0 
LD=0 
ISGN=1 

IF(L3-35)317,317,323 
323 LI=LI+1 
L2 = 7 

L3=L3-36 

G0TO(31,26,45,47),L3 

317 IF(IABS(IS)-1)99,338,324 

338 15=1 
L3=L3-33 

GOTO(320,337) , L3 

320 IS=-2 
L2 =3 

337 M9=M9+1 

IF(M9-M) 339,339, 326 

339 GOTO(47,311) ,L3 

321 IF(M0DMM(M9) /1 00-4) 34,35,34 

45 L2=L2+1 

46 L2=L2+1 

47 CONTINUE 

IF(Ll-6 >330,336,99 

336 LI=LI1 
IND=4 

IF( IABS( IS >-1)99,60,324 
330 IF( ( 10*Ll+L2)/ll-3>98,48,49 

48 INC(2,M)=0 

49 1=1 

IF! Ll-5) 50,7 ,99 

50 Z(M)=FLOAT( IZ) 

Z(M)=Z( M)/10.0**= (LD) 

GOTO 6 

51 N1=N+1 

IFIMTRANd)) 99,95,202 

202 IF( I\1LST( J + 8)-I SYMB(26) ) 204,52,204 
204 IF( INLST(J+8>- ISYMB(13 ) >55,201, 55 
201 DO 203 J=l,2 

DO 200 1=1,19 
K=IRAY( I, J ) 

200 IRAY( I ,J)=ISYMB(K) 

203 WRITE ( 6, 116){ IRAY( I ,J ) , 1=1, 19) 

52 IPRIN=1 

55 NA=0 

DO 58 1=1, N,1 
IF(NODE( I)-l ) 58,56,57 

56 WRITE (6,112) I 

57 NA=NA+1 
NODE(NA)=I 

58 CONTINUE 
IF(NASA-2)70,71,70 

70 NASA=4 
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f 



71 

59 



91 

60 

95 

96 
90 

97 
72 

98 

99 

326 

328 

324 

100 

101 

106 

107 

108 

109 

110 



GO TO 66 

NASA=3 

IND=2 

IF(N-NA) 99,60,59 
WRITE (6,111) (NODEd) 
WRITE (6,110) 

GO TO 60 

IND=7 

RETURN 

WRITE (6,115) 

GO TO 90 

114) 



1=1, NA, 1) 



MOVFL 



106) 



INLST 



J 

INLST 



111 



NAS AP 
CONTINUE 



PROBLEM 
ON THE 



112 

] 

1 14 

115 
116 

117 

119 

304 

1 

325 

327 

] 

329 

] 

376 

998 



1 



WRITE (6 
NASA=4 
GOTO 66 
WRITE (6 
WRITE (6,116) 

NASA=5 
GO TO 49 
WRITE (6,108) 

WRITE (6,116) 

IF(Ll-6)72,66, 99 
WRITE (6,107 ) 

GO TO 90 
WRITE (6,327) 

GOTO 66 
WRITE (6,329) 

WRITE (6,325) 

GOTO 66 

FORMAT ( Ell. 3, 9X, El 1.4) 

FORMAT (80A1 ) 

F0RMAT(24H FUNCTIONS NOT AVAILABLE) 

FORMAT! 27H MACHINE LANGUAGE BREAKDOWN) 

F0RMAT(32H INPUT CODING ERROR IN COLUMN, 13) 
FORMAT! 39H USE HEADING LABELED 
FORMAT!/, IX, 39HTHE ANALYSIS WILL 
,5HPTI ON, 

26H THAT THE LIST IS CORRECT. 

,9HPING NODE, 

39H NUMBERS IS NOT A RECOMMENDED PRACTICE.//) 
FORMAT! IX, 41HTHE FOLLOWING NODE NUMBERS ARE MENTI 
,6HIN THE, 

12H INPUT LIST , / , IX , I 2 , 19 ( IH , , I 2 , IX ) ) 

FORMAT!//, IX, 5HN0DE ,I2,27H HAS ONLY ONE ELEMENT 
IIHCTED TO IT. ,// ) 

FORMAT! IX, 12HN0 MORE THAN,I3,18H ELEMENTS ALLOWED 
FORMAT ( IX, 24HTHERE IS NO OUTPUT LIST.) 

FORMAT (IX, 80A1 ) 

FORMAT! IX, 1116 ) 

FORMAT (1 HI ) 

FORMAT! //IX, 39HM0D1 FY AND SYMBOLIC SOLUTION MUST 
,9HW EXECUTE//) 

F0RMAT(//1X,39HF0R LOOPS YOU MUST SPECIFY — (START 
, 4H I NC ) ) 

F0RMAT(//1X,39HT00 MANY ELEMENTS OR VALUES IN MOD 
,8HTATEMENT) 

FORMAT! //IX, 39HSTART GREATER THAN END AND/OR INC 
,4HARGE) 

FORMAT (/1X,27HTHE TRANSFER FUNCTION INPUT/ 
1X,33HVARIABLE IS A DEPENDENT QUANTITY.) 

STOP 

END 



ASSUM 

/ , 1X,13HH0WEVER, SKIP 

ONED 

CONNE 

.) 

FOLLO 
,END, 
IFY S 
TOO L 



FUNCTION ITRAN(I) 
ITRAN=1+I/16777216 
IF! ITRAN)1,2,3 

1 ITRAN=I TRAN+63 
RETURN 

3 ITRAN= ITRAN-31 
RETURN 

2 WRITE(6,4) 
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" V#=- 



10 ^ 



I 



I 



4 



RETURN 

FORMAT! //IX, 15HMACHINE FAILURE) 

STOP 

END 



SUBROUTINE NAS4 

INTEGER SENSFL,FLAG,EXPNS,SIGN,WHERE,EXPDS,S IGD, 
IPARAM, PARTLN, PARTED 

DIME NS ION JZ(5) , AZ( 5) , LNUM ( 23, 3 ) , L DEN ( 23 , 3 ) , ILGTH( 23) , 
IIFMTK 18) ,IFMT2( 12) 

COMMON/SEN/PARAM! 3) , SEN SFL , FRMAX , FRMI N 
COMMON IND,NASA,MTRN1, MTRN2 ,MTRN3,HD, IJTAG(15,2),I 
1PRIN,Z(42) ,LL, M0DNM(23) , VALMOD ( 23 , 10) ,IL00P(23),LI,LIC 
2,M9, ZZ(23) ,NAM(23,3 ), ITAG( 5,8) ,NMBR(6) , 

3M,SNUM( 3 0,5) ,SDE N( 3 0 ) , I OD ( 30 ) , I 0 ( 5 ) , NUM ( 5 ) 

5, XG( 5) , IP( 30,5 ) , ID,NUMZ,DMY( 23) 

4, NAME (23, 2 ) , M\/ 1 , 1ST ( 40 ) , L IM ( 6 ) , I TWO (23) ,LOOP( 10000) 

6, WHERE,PARTLN( 200,23) ,EXPNS 

7 (200),NFN( 200) , PARTED! 200, 23 ) , E XPD S( 200) , 

8SIGD(200) ,NFD( 200) , LPN,LPD, S IGN( 200 ) 

DATA IL3TH/1H1, 1H2, 1H3 , 1H4 , 1H5 , 1H6 , 1H7 , 1H8 , 1H9 , 2H10 , 
12H11,2H12,2H13,2H14,2H15, 2H16, 2H17, 2H18, 2H19, 2H20,2H21 
2,2H22,2H23/ 

DATA IFMT1/4H! 1 X , , 4HI 5 , 1 ,4H8 , 3X ,4H , A 1 , , 4H 1 H ( , , IHl , 
14H(A1, ,1H2,4HI1),,4H3H /,2H , , IHi, 4H( A 1 , , 4H2I 1 ) , 
24H,6H),4H S * , 4H^ , I 3 , IH ) / 

DATA IFMT2/4H! IX , , 4HI 5 , I , 4H8 , 3X , 4H , A 1 , , 4H 1H( ,,1H1, 
14H(A1, ,1H2,4HI 1) , ,4H6H) ,4HS **,4H,I3)/ 

DATA IM/»- '/,IB/' '/,IA0C/'A1, '/,ITW/'2 '/ 

MDD=MTRN1+1 

IF(IPRIN-1)1,1,23 

1 IERR=0 
IF(LIC)60, 2, 5 

2 DO 3 N1=1,MVI 

3 ZZ(N1)=Z(N1) 

IFdLOOPd )-10 00 )5,16,60 

16 DO 19 Nl=l ,MTRN1 
DO 19 L = l, 2 

19 IJTAG! N1,L)=IJTAG(N1,L )-1000* ( IJTAGI N 1, L ) / 1000) 

IF( IPRIN)60,13,23 
23 1 = 0 

LPN=0 

DO 1024 K2=1,MDD 
N7=LIM( K2) 

WHERE=0 

IF(K2-MDD) 50,51,60 

50 WRITE (6,135) ( ITAG ( K2 , K ) , K= 1 , 8 ) 

WRITE (6,160) 

J=-l 

IF(NMBR(K2 ) )601,52,600 

600 JS=IM 
GOTO 602 

601 JS=IB 

602 WRITE (6,162) JS 
GOTO 52 

51 WRITE (6,138) 

WRITE (6,160) 

JS = IB 

WRITE (6,162) JS 
SIGD( 1)=JS 
NFD( 1)=1 
LPD=1 
J = 1 

52 N8=l 
IF(N7)60,55,55 

55 WRITE (6,100) 

WRITE (6,141) 

GO TO 1024 

56 DO 24 N9=1,N7 
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FLAG=0 
KK=0 
KKK=0 
ISL0P=0 
JS = I B 
1 = 1+1 

LK1 = L00P ( I ) 

LK2=IABS(LK1 ) 

IF(J#LK1)44,60,25 
44 JS=IM 

25 N1=N8 

26 IF(LK2-ITW0(N1) ) 27,28,28 

27 IF(N1-N8 )60,3D ,29 

30 N8=N8+1 
GOTO 25 

28 ISLOP=ISLOP+IST (N1 ) 

NM=NAME(N1,1)/10D 

IF( (NM + 3)=?(NM-4)*{NM-5) ) 32, 31, 32 

31 KKK=KKK+1 

DO 230 K=l,3 

230 LDEM(KKK,K ) = NAM(N1,K) 

GO TO 33 

32 KK=KK+1 

DO 231 K=l,3 

231 LNUM(KK,K) =NAM (N1 ,K) 

33 LK2=LK2-ITW0(\| 1) 

29 N1=N1+1 

IF (Nl-MVI-1) 26 ,499,60 

499 IF(SENSFL-l) 34,500,60 

500 JJ=1 

502 IF( K2-MDD) 509,799,60 

509 LPN=LPN+1 
IF(LPN.GT.200) WRITE(6,164) 

IF(KK.EO.O) GO TO 531 

510 JJJ=1 

DO 511 K=l,3 
IF(K.GT.l) JJJ=4 

CALL PACK(PARTLN(LPN,J J) ,K,LNUM( JJ,K) ,JJJ) 

511 CONTINUE 

IF(LNUM( JJ, 1) .EO.PARAMC D.AND. LNUM( JJ,2) . E Q. PAR AM (2) .A 
1ND.LNUM(JJ,3).E0.PARAM(3)) GO TO 519 
CALL PACK( PARTLN(LPN,JJ) ,4,0,4) 

GO TO 530 
519 FLAG=1 

CALL PACK(PARTLN(LPN,JJ) ,4,3,4) 

530 JJ=JJ+1 
IF(JJ-KK-1)510, 531,60 

531 IF(KKK.EQ.O) GO TO 554 
LX=0 

JJ=JJ-1 
534 JJ=JJ+1 
LX=LX+1 
JJJ = 1 

DO 532 K=l,3 
IF(K.GT.l) JJJ=4 

CALL PACK(PARTLN(LPN, JJ) ,K,LDEN(LX,K), JJJ ) 

532 CONTINUE 

533 IF( LDEN( LX,1 ), EQ.PARAMI 1) .AND.LDEN(LX,2) .EQ.PARAM(2) .A 
IND. LDEN( LX,3) . EQ. PARAM(3) ) GO TO 540 

CALL PACK(PARTLN(LPN, JJ ) ,4,6,4) 

GO TO 542 
540 FLAG=1 

CALL PACK! PARTLN(LPN,JJ) ,4,7,4) 

542 IFILX-KKK) 534,555,60 

554 JJ=JJ-1 

555 WHERE=1 
GO TO 34 

799 LPD=LPD+1 
IF(LPD.GT.200 ) WRITE(6,164) 

IF(KK.EO.O) GO TO 831 

800 JJJ=1 

DO 812 K=l,3 
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IF(K.GT.l) JJJ=-Jf 

CALL PACK(PARTLD(LPD, J J) , K , L NUM ( J J , K ) ,JJJ) 

812 CONTINUE 

IF(LNUM{ JJtD.EO.PARAMd) . AND. LNUM( J J ,2 ) .EQ. PAR AM (2) .A 
IND.LNUMI JJ, 3) .EO. PARAMO) ) GO TO 819 
CALL PACK(PART LD(LPD,JJ ) ,4,0,4 ) 

GO TO 830 
819 FLAG=1 

CALL PACK{PARTLD(LPD,JJ) ,4,3,4) 

830 JJ=JJ+1 

IF{ JJ-KK-1 >800,831,60 

831 IFIKKK.EO.O) GO TO 854 
LX=0 

JJ=JJ^1 
834 JJ=JJ+1 
LX=LX+1 
DO 832 K=1 ,3 
JJJ=1 

IF(K.GT.l) JJJ=4 

CALL PACK(PARTLD(LPD, JJ) ,K,LDEN(LX,K), JJJ ) 

832 CONTINUE 

83 3 IF{ LDEN( LX, 1 ). EO.PARAMI 1 ) . AND. LDEN { L X, 2 ) . EQ. P ARAM ( 2 ) . A 
1ND.LDEN( LX,3). EO.PARAMO) ) GO TO 840 
CALL PACK(PARTLD(LPD, JJ ) ,4,6,4) 

GO TO 842 
840 FLAG=1 

CALL PACKI PARTLDI I , JJ ) ,4,7,4) 

842 IF(LX-KKK) 834,855,60 

854 JJ=JJ-1 

855 WHERE=2 

34 IF (KK) 60,42,43 

42 KK=1 
LNUM{1,1 )=IB 
LNUMd ,2) =18 
LNUM{ 1, 3)= 1 
IFMTl (8 )=1A0C 

43 IF{ KKK) 60,40,4 1 

40 IFMT2(6)=ILGTH(KK) 

WRITE(6,I FMT2) I ,LK1 , JS , ( ( LNUM ( K , K1 ) , Kl= 1 , 3 ) , K= 1 , KK ) , IS 
ILOP 

GO TO 901 

41 IFMT1(6)=ILGTH(KK) 

IFMTK 12) = ILGTH{ KKK) 

WRITE! 6, IFMTl )I,LK1,JS, ( (LNUM{K,K1) ,K1=1 , 3) , K=1 ,KK) , ( ( 
1LDEN(K, Kl) ,K1 = 1 ,3) ,K=1 ,KKK) , I SLOP 
IF{ I FMTK 8 )-ITW) 232,901,232 
232 IFMT1(8)=ITW 

901 IF{ WHERE-1 >24,902,904 

902 EXPNS! LPN)=ISL0P 
SIGN(LPN)=JS 
NFN( LPN) =J J 
IF(FLAG)60,905,932 

904 EXPOS! LPD)=ISLOP 
SIGD! LPD) =JS 
NFD!LPD)=JJ 
IF!FLAG)60,903 ,932 

905 DO 906 K=1,JJ 
ITEST=0 

CALL PACK! IT ES T , 4, P ARTLN ! LPN , K ) , 4 ) 

IF! ITEST-5)960,60,961 

960 CALL PACK!PARTLN!LPN,K) ,4,4,4) 

GO TO 906 

961 CALL PACK! PARTLN! LPN, K) ,4,9,4) 

906 CONTINUE 
GO TO 932 

908 DO 931 K=1,JJ 
ITEST=0 

CALL PACK! I TEST, 4, PARTL D! LPD, K ) ,4) 
IF!ITEST-5)962,60,963 

962 CALL PACK!PARTLD!LPD,K) ,4,4,4) 

GO TO 931 

963 CALL PACK!PARTLD!LPD,K) ,4,9,4) 
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931 CONTINUE 

932 IF(N9.E0.N7) 30 TO 909 
GO TO 24 

909 IF(WHERE-1)24,930»930 

930 CALL SENS 
24 CONTINUE 
1024 CONTINUE 

IPRIN=IPRIN-1 

IF( IPRIN)13,13,89 ' 

5 LIC=LIC+1 
WRITE (6,163) 

DO 12 N1=1,M9 

IF( LIC-IABS( I LOOP! Nl) ) )4,4,12 
4 DO 11 N7=1,MVI 

IF{ MODNMINl )-I ABS(NAME(N7,1 ) ) )11,6, 11 

6 IF( ILOCP(Nl) 1 8,60,7 

7 ZZ(N7) = VALM0D{ N1,LIC) 

GOTO 9 

8 ZZ(N7) =VALMOD( N 1 , 1 ) +FL OAT { L I C-1 ) »!^ VALMOD ( Nl ,3) 

9 NM=NAME(N7, 1 )/ 100 

WRITE (6,129) (NAM(N7,K) ,K=1,3),ZZ(N7) 

IF( (NM + 3)'!'(NM-4)*(NM-5) ) 12, 10,12 

10 ZZ(N7)=1./ZZ(M7) 

GOTO 12 

11 CONTINUE 
WRITE (6,130) 

I ERR=1 

MVI 1=MVI+1 
DO 15 N7=MVI1, M 

IF( M0DNM(N1)-I ABS( NAME(N7,1) )) 15, 14, 15 

14 WRITE (6,133) MODNM(Nl) 

GOTO 12 

15 CONTINUE 

WRITE (6,132) MODNM(Nl) 

12 CONTINUE 

IF( IERR)60,13, 66 

13 IF( ILOOP (1 )-lD00)17, 18 ,60 
18 WRITE (6,101) 

17 DO 400 N 1=1, 30 
SDEN(N1 ) =0 . 

DO 400 K2=1,MTRN1 

400 SNUM(N1,K2)=0. 

DO 603 K2=1,MTRN1 
603 SNUM( 16,K2)=N^BR(K2) 

SDEN( 16)=1. 

I =0 

DO 1412 K2=1,MDD 
N8 = l 

N7=LIM(K2) 

IF(N7)60,1412, 397 
397 DO 412 N9=1,N7 
1 = 1 + 1 
ISL0P=0 
VL00P=1. 

LK1 = L0GP(I ) 

LK2= IABS(LK1 ) 

SGN=LK1/LK2 

401 N1=N8 

402 IF(LK2-ITW0(N1 ) ) 403,405,405 

403 IF (Nl-N8)60,404,406 

404 N8=N8+1 
GOTO 401 

405 ISLOP=I SL0?+IST(N1) 

VL00P=VL00P*ZZ(N1) 

LK2=LK2-ITW0(N1 ) 

406 N1=N1+1 
IF(N1-MVI-1)402,407,60 

407 JS=16-ISL0P 
IF(MDD-K2) 60,410,408 

40 8 SNUM( JS ,K2) = SNUM( JS,K2)+VL00P=«-'SGN 
GOTO 412 

410 SOEN( JS)=SDEN( JS)+VLOOP*SGN 
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412 CONTINUE 
1412 CONTINUE 

DO 69 K=1,MTRN1 

JZ(K)=1 

AZ(K) =0. 

DO 68 J=l, 30 
IF(SNUM( J,K) )67, 68,67 
67 IF ( AZ(K) ) 201 »202 ,201 
202 AZ(K)=-SNUM( J,K) 

201 JZ(K)=J 

68 CONTINUE 

69 CONTINUE 
SIGNC=0. 

JI=0 

DO 71 J=l,30 
IF(SDEN( J) )97, 71,97 

97 IF( SIGNC) 70, 98, 70 

98 SIGNC=SDEN(J ) 

70 JI=J 

71 CONTINUE 
J = 0 

DO 73 K=1,MTRN1 
IF(JZ(K)-J) 73,73,72 

72 J=JZ(K) 

73 CONTINUE 
IF(JI-J) 75,75,74 

74 J=JI 

75 DO 83 K=1,MTRN1 
XG(K)=0. 

IF ISIGNC)90,91 ,90 
90 XG( K) =AZ (K) /SIGNC 

91 WRITE (6,136) XG(K) 

WRITE (6,135) ( ITAG(K, I) ,1=1,8) 

NU = 0 
I0(K)=0 
DO 82 L=1,J 
ANUM=NU 

IF(ABS(SNUM(L,K) )* ( ANUM-. 5 ) ) 80 , 82 , 81 

80 IO(K)=J-L 

81 NU=NU+1 
IP(NU,K )=J-L 

SNUM(NU.K) =-SNUM( L, K) /AZ(K) 

WRITE (6,137) SNUM( NU,K) ,I P(NU,K) 

82 CONTINUE 
NUM( K) =NU 

IF(NUM(K)) 60,92,83 

92 NUM(K)=1 

WRITE (6,137) SNUM(1,K) ,IP(1,K) 
WRITE (6,141) 

83 CONTINUE 
WRITE (6,138) 

NUMZ=0 

ID=0 

DO 86 K=1,J 
ANUM=NUMZ 

IF(ABS(SDEN(K) )* (ANUM-.5) ) 84, 86, 85 

84 ID=J-K 

85 NUMZ=NUMZ+1 
IOD( NUMZ)=J-K 

SDEN( NUMZ) =SDEN(K) /SIGNC 

WRITE (6,137) SDEN(NUMZ),IOD(NUMZ) 

86 CONTINUE 
IF(NUMZ) 60,94,87 

94 NUMZ=1 

WRITE (6,137) SDEN(l) ,I0D(1) 

WRITE (6,141) 

87 IF(LI-LIC)150, 150, 152 
150 LIC=0 

152 IF(MTRN1-MTRN2 )96, 96, 95 

95 IF(MTRN1-MTRN3 ) 88,88, 153 

153 IF(LIC) 60,89,5 

96 IND=5 
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GO TO 900 

88 IND=6 

GO TO 900 

60 WRITE (6,134) 

66 LIC=0 

89 IND=1 
900 RETURN 

100 FORMAT! //IX, 16HIDENTICALLY ZERO) 

101 FORMAT ( ///IX, 1 6HN0M INAL SOLUTION) 

129 FORMAT! /IX, lOHMODI FYING ,A1,2I1,4HT0 ,E11.4) 

130 FORMAT! //IX, 26HERR0R IN ATTEMPT TO MODIFY) 

132 FORMAT ! IX, 12HBRANCH CODED, 15, 22H — NOT IW ORIGINAL NETW 
1,3HDRK) 

133 FORMAT ! IX, 12HBRANCH CODED, 15, 22H — CANNOT MODIFY INDEPE 
1 ,13HNDENT SOURCES) 

134 FORMAT! //IX, 24HMACHINE LANGUAGE FAILURE! 

135 FORMAT ! /IX, 24HNUMERAT0R POLYNOMIAL OF , 2A 1 , 21 1 , IH / , 2A 1 

1 , 211 ) 

136 FORMAT! //IX, 4HK = ,E11.4) 

137 FORMAT! IX, Ell. 4, lOH TIMES S>?==«',13) 

138 FORMAT! //IX, 22HDEN3MI NATOR POLYNOMIAL) 

141 FORMAT !/, IX, 31HCHECK YOUR CIRCUIT CODING LIST., 
+/,1X,26HTHIS ANSWER IS DEGENERATE.) 

160 FORMAT! /3X,28HL0C CODE LOOP PRODUCT) 

162 FORMAT ! 17X,A1, 13H! 1 ) S 0) 

163 FORMAT!/) 

164 FORMAT! • «,• WARNING: NUMBER OF UNCANCELLED LOOP PRODU* 

1, 'CTS EXCEEDS STORAGE AVAILABLE !200). »,/• •,'YOU MU* 

2, ' ST CHANGE DIMENSION STATEMENTS IN NAS4 AND SENS SUB* 

3, 'ROUTINES FOR VARIABLES - ',//* • , • PARTLN , EXPNS , NFN* 

4, • , PARTED, EXPDS, SIGD,NFD, S IGN, RNUMP,RNUMS, RDENS,RDENP* 
5) 

STOP 

END 

C PACKING SUBROUTINE 



SUBROUTINE PACK! A, I ,B , J) 
LOGICAL*! A! 1) ,B( 1) 
A!I)=B! J) 

RETURN 

END 



SUBROUTINE SENS 

COMPLEXES CFRE0,CTEMP,CTEMP1 ,CTEMPB 
INTEGER EXPNS, SIGN , EXPDS, SI GD, WHERE ,PARAM, PARTED, 
IPARTLN 

DIMENSION RNUMP!200) ,RNUMS!200) ,RDENS!200) ,RDENP!200), 
1 ITEMP! 23 ) , lEXPN! 50 ) , IE XPD! 50) , PARTI! 501 ,PART2! 50) 
C0MM0N/SEN/PARAM!3 ) , S ENS FL , FRM AX , FRM IN 
COMMON IND,NASA,MTRN1,MTRN2,MTRN3,MD,IJTAG!15 ,2) , I 
IPRIN, Z!42) , LL, M0DNM!23 ) , V ALMOD! 23 , 10) , I LOOP ! 23) , L I , L I C 
2,M9,ZZ!23) , NAM ! 23 , 3 ) , I TAG ! 5 , 8 ) ,NMBR!6), 

3M,SNUM! 3 0, 5) ,SDEN!30) , lOD! 30) , I 0! 5) ,NUM!5) 

5 , XG(5 ) , IP!30,5 ) , ID,NUMZ, DMY ( 23 ) 

4, NAME! 23,2) ,MVI ,IST!40) ,LIM!6) , I TWO (23) , LOOP (10000) 

6, WHERE, PARTLN! 200, 23 ) , E XPN S 

7(200) , NFN! 200) , P ARTL D! 200 , 23 ) , EXPDS (200) , 

8SIGD! 200) , NFD! 2 00) ,LPN , LPD , SI GN ! 2 00 ) 

DATA INEG/*- '/,IBLNK/' ’/ , PARTl/50*0. 0/, PART2/50 
1*0.0/ 

DO 2 K=1,MVI 
1TEMP!K)=NAM!K,1) 

DO 2 L=2,4 

CALL PACK!ITEMP!K) ,L,0,4) 

2 CONTINUE 

IF! WHERE-1)60,200,19 
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19 CALL PACK(PARTLD{1,1) ,3,1,4) 

CALL PACK{PARTLD( 1 , 1) ,4,4,4) 

NFD( 1) =1 
GO TO 400 

200 DO 320 I =1 ,LPN 
RNUMP( I ) =1 
RNUMS( I )=1 
KP=NFN( I ) 

DO 320 J=1,KP 

IC0MP1=0 

ICOMP2=0 

IC0MP3=0 

ICOMP4=0 

CALL PACK ( I COM PI ,1 , PART LN ( I , J ) , 1 ) 

CALL PACK( IC0MP2,4,PARTLN( I , J) ,2) 

CALL PACK( IC0MP3,4,PARTLN( I,J),3) 

KL = 0 

205 KL=KL+1 

IF{ ICOMPl.EO.ITEMP(KL) ) GO TO 210 
GO TO 225 

210 IF( IC0MP2.E0.NAM(KL,2) ) GO TO 220 
GO TO 225 

220 IF( IC0MP3.E0.NAM( KL,3) ) GO TO 240 

225 IF(KL.LT .MVI ) GO TO 205 

240 CALL PACK(IC0MP4,4,PARTLN(I,J),4) 

ICOMP4=ICOMP4+1 

GO TO (250,60, 60,285,280,60, 290, 300,60, 310) , IC0MP4 
250 RNUMP( I )=RNUMP( I )/Z(KL) 

GO TO 285 
280 RNUMP(I)=0 
285 RNUMS( I ) =RNUMS ( I ) /Z(KL) 

GO TO 316 

290 RNU.MPd )=RNUMP(I)*Z(KL) 

GO TO 315 

300 RNUMP( I )=RNUMP { I )=!=Z{KL )*=^=2 
RNUMPC I ) =-RNU,MP( I ) 

GO TO 315 
310 RNUMP(I)=0 

315 RNUMS( I ) =RNUMS ( I )* Z(KL) 

316 IF(SIG.M( D.EO.INEG) GO TO 317 
GO TO 320 

317 RNUMP( I ) =-RNUMP( I ) 

RNUMS( I )=-RNUMS( I ) 

320 CONTINUE 
GO TO 30 

400 DO 2520 1=1, LPD 
KP = NFD( I ) 

ROENP{ I ) =1 
RDENS( I )=1 

IF{ I.NE.l) GO TO 402 
RDENP( I )=0 
GO TO 2520 
402 DO 520 J=1,KP 
ICOMP1=0 
IC0MP2=0 
IC0MP3=0 
IC0MP4=0 

CALL PACK{ IC0MP1,1,PARTLD( I, J), 1) 

CALL PACK(IC0MP2,4,PARTLD(I , J) ,2) 

CALL PACK( IC0MP3,4,PARTLD( I ,J) ,3) 

KL=0 

405 KL=KL+1 

IF{ ICOMPl.EO.ITEMP(KL) ) GO TO 410 
GO TO 425 

410 IF( IC0MP2.E0.NAM(KL,2) ) GO TO 420 
GO TO 425 

420 IF( IC0MP3.E0.NAM(KL,3) ) GO TO 440 

425 IF(KL.LT.MVI ) GO TO 405 

440 CALL PACK! IC0MP4,4,PARTLD( I, J) ,4) 

IC0MP4=IC0MP4f 1 

GO TO (450,60,60,485,480,60,490,500,60,510) ,IC0MP4 
450 RDENP( I )=RDENP ( I )/Z(KL ) 
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GO TO 485 
480 RDENP(I)=0.0 
485 RDENS ( I )=RDENS ( I )/ Z(KL ) 

GO TO 516 

490 RDENP( n=RDENP ( I )’^=Z(KL) 

GO TO 515 

500 RDENP( I )=RDENP ( I )*Z(KL)**2 
RDENP ( I )=-RDENP( I ) 

GO TO 515 
510 RDENP(I)=0.0 

515 RDENS( I )=RDENS ( I )*Z(KL ) 

516 IF{ SIGD( I) .EO. INEG) GO TO 517 
GO TO 520 

517 RDENP( I ) =-RDENP(I ) 

RDENS( I ) =-RDENS( I ) 

520 CONTINUE 
2520 CONTINUE 
J = 0 

DO 610 K=1,LPD 
DO 610 I =1 ,LPN 
J = J + 1 

IEXPN( J)=EXPNS ( I )+EXPDS(K) 

IF(ABS( RDENS (K)*RNUMP( I)-RNUMS ( I )*RDENP(K) ) .LE.l.OE-75 
1) GO TO 1520 
GO TO 1521 

1520 J=J-1 

GO TO 610 

1521 IFU.EO.l) GO TO 1620 
JHIN 1=J-1 

DO 1610 KL=ltJMINl 

IF(IEXPNU).EQ.IEXPN(KL) ) GO TO 608 
GO TO 1610 

608 PARTI (KL)=PART1 (KL ) +RDENS ( K )=«^RNUMP( I )-RNUMS( I )*RDENP(K 
1) 

J = J-1 
GO TO 610 
1610 CONTINUE 

16 20 PARTK J )=RDENS( K)*RNUMP( I )-RNUMS( I I-i'RDENPI K) 

610 CONTINUE 

IF(J.LE.l) GO TO 640 
JJ = 1 

631 K=JJ+1 

IF(ABS( PARTK J J) ) . LE. 1. OE-75) GO TO 636 

632 IF( lEXPN(JJ) .GT.IEXPN(K) ) GO TO 634 
TEMP1 = PART1 ( JJ ) 

TEMP2=I EXPN( JJ ) 

PART1(JJ)=PART1(K) 

IEXPN( JJ) = IEXPN( K) 

PARTKK )=TEMP1 
IEXPN(K)=TEMP2 
634 K=K+1 

IF(K.GT.J) GO TO 636 
GO TO 632 

636 IF( JJ.EO«( J-1) } GO TO 640 
JJ=JJ+1 

638 GO TO 631 

640 WRITE(6,900) ( I TAG ( 1 , K ) , K=1 , 8 ) , ( P AR AM( K ) » K=1 , 3 ) 

IF(J 160,641,642 

641 WRITE(6,906) 

GO TO 643 

642 SET=0 

DO 710 K=1,J 

IF(ABS(PART1(K) ) .LT.l. OE-75) GO TO 710 
SET=1 

700 WRITE(6,901) P ARTl ( K ) , I EXPN ( K ) 

710 CONTINUE 

IF(SET.EO.O) /JRITE(6,906) 

643 J1=J 
J = 0 

DO 611 K=1,LPD 
DO 611 I =1 ,LPD 
J = J + 1 
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IEXPD( J )=EXPDS ( I )+EXPDS(K) 

IF(ABS(RDENS(K)=<'RDENS( I)).LE.1.0E-75) GOTO 1711 
GO TO 1712 

1711 J=J-1 

GO TO 611 

1712 IF(J.EO.l) GO TO 1612 
JMIN1=J-1 

DO 1611 KL=1,JMIN1 

IF( IEXPD( J ) .EQ. lEXPD(KL) ) GO TO 609 
GO TO 1611 

609 PART2{ KL ) =PART2( KL)+RDENS( K)*RDENS ( I ) 

J = J-1 
GO TO 611 

1611 CONTINUE 

1612 PART2( J)=RDENS(K)*RDENS( I ) 

611 CONTINUE 

IF(J-1)60,749,740 

740 JJ=1 

741 K=JJ+1 

IF(ABS(PART2(J J) ).LE.1.0E-75) GO TO 746 

742 IF(IEXPD(JJ).GT.IEXPD(K) ) GO TO 744 
TEMP1 = PART2( JJ ) 

TFMP2=I EXPD( JJ ) 

PART2( JJ)=PART2(K) 

IEXPD( JJ)=IEXPD(K) 

PART2(K )=TEMP1 
lEXPD(K) =TEMP2 
744 K=K+1 

IF(K.GT.J) GO TO 746 
GO TO 742 

746 IF(JJ.EO.{ J-1) ) GO TO 749 
JJ=JJ+1 

748 GO TO 741 

749 WRITE(6»902) 

DO 750 K=1,J 

IF(ABS(PART2(K) ) .LT.l.OE-75) GO TO 750 
WRITE(6,901) PART2(K), lEXPD(K) 

750 CONTINUE 

IF(FRMAX.EO.O.O.AND.FRMIN.EQ.O.O) GO TO 30 
FRINC= (FRMAX-FRM IN )/ 199.0 
FREO=FRMIN 
WRITE! 6, 904) 

DO 763 KK=1,200 

FREO 1=2. 0=^=3. 14 15 92 653 5=)' FREQ 

CFRE0=CMPLX(0.0, FREQl) 

CTEMP=(0.0,0.0 ) 

DO 755 K=1»J1 

CTEMP1=CHPLX(PART1(K) , 0.0) 

IF (lEXPN(K) .LT.O) GO TO 754 
CTEMP = CTEMP+CTEMP1>!'CFREQ**I EXPN( K) 

GO TO 755 

754 CTEMP = CTEMP + CTEMP1/CFREQ=«'*(IABS(IEXPN(K) ) ) 

755 CONTINUE 

CTEMPB=(0.0,0.0 ) 

DO 760 K=1,J 

CTEMP1=CMPLX(PART2( K) , 0. 0) 

IFdEXPD(K).LT.O) GO TO 759 
CTEMPB=CTEMPB+CTEMP1*CFREQ=^'*I EXPD(K) . 

GO TO 760 

7 59 CTEMPB=CTEMPB+CTEMP1/CFREQ**(I ABS( IEXPD(K ) ) ) 

760 CONTINUE 

IF(CABS(CTEMPB ).LT.1.0E-74) GO TO 1762 

CTEMP=CTEMP/CTEMPB 

RMAGN=CABS(CTEMP ) 

ANGL E= AT AN2 (AI MAGI CT EMP }» REAL! CTEMP) ) 

762 WRITE(6»905) F RE Q , RMAG N » ANGLE 
1762 FREO=FREO+FRINC 

763 CONTINUE 
GO TO 30 

900 FORMAT! • ' , / IX » 3 IHNUMERATOR POLYNOMIAL : SENSITIV 

1 , ' ITY OF • »2A1 ,211 » 

21H/,2A1,2I1,2X,4HT0 ,A1,2I1) 
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901 F0RMAT(« ' , / , IX , El 1 .4, lOH TIMES 

902 FORMATC ‘ ,////, IX , 27HDEN0M IN ATOR POLYNOMIAL : S 

1» 'ENSITIVI TY* ) 

903 FORMATC ‘ MACH INE FAILURE') 

904 FORMATC ',////• ' , • EV A LUAT I ON OF SENSITIVITY', 

1///,' ', 'FREQUENCY' 

2 ,10X, ' MAGNITUDE' , lOX, ' PHASE! IN RADIANS)',/,' ') 

905 FORMATC ' ,E 1 1 . 4 , 8X , El 1 . 4 , 8X , El 1 . 4 ) 

906 FORMAT!' ',//,' ' , ' I DENT ICALLY ZERO',//,' ») 

60 WRITE{6,903) 

30 RETURN 
END 
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